Non-universal size dependence of the free energy of confined systems near criticality 
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The singular part of the finite-size free energy density fs of the 0(n) symmetric ip"^ field theory in 
the large-n limit is calculated at finite cutoff for confined geometries of linear size L with periodic 
boundary conditions in 2 < d < 4 dimensions. We find that a sharp cutoff A causes a non-universal 
leading size dependence fs ~ h^~^L~^ near Tc which dominates the universal scaling term ~ L^"*. 
This implies a non-universal critical Casimir effect at and a leading non-scaling term ~ L^^ of 
the finite-size specific heat above Tc. 
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The concept of universal finite-size scaling has played 
an important role in the investigation of finite-size effects 
near critical points over the last decades . Consider 
the free-energy density f{t,L) of a finite system at the 
reduced temperature t — {T—Tc)/Tc and at vanishing ex- 
ternal field in a d-dimensional cubic geometry of volume 
L'^ with periodic boundary conditions. It is well known 
that, for small t, the bulk free energy density = f{t, oo) 
can be decomposed as 



hit) = fbsit) + fo{t) 



(1) 



where fts — A^\t\'^~°' denotes the singular part of /f, 
above (+) and below (— ) Tc, apart from bulk corrections 
to scaling, and where the regular part fo{t) of fb can be 
identified unambiguously. 

According to Privman and Fisher ^,|^ the singular 
part of the finite-size free-energy density may be defined 

by 



fs{t,L) = f{t,L)-fo{t) 



(2) 



where /o is independent of L. The finite-size scaling 
hypothesis asserts that, below the upper critical dimen- 
sion d = A and in the absence of long-range interactions, 
fs{t,L) has the structure 



(3) 



where T{x) is a universal scaling function and ^ = 
is the bulk correlation length. Both ^ and L are 
assumed to be sufficiently large compared to microscopic 
lengths (for example, the lattice spacing a of lattice mod- 
els, the inverse cutoff A^^ of field theories, or the length 
scale of subleading long-range interactions). Eqs.(|l|) - 
are expected to remain valid also for non-cubic ge- 
ometries where the scaling function J-{x) depends on the 
geometry and on the universality class of the bulk criti- 
cal point but not on a or A and not on other interaction 
details |-§. 

As a consequence, universal finite-size scaling proper- 
ties are generally believed to hold for quantities derived 



s, such as the critical Casimir force in film geome- 
10 1 and the finite-size specific heat Cs, in apparent 



from 
try 

agreement with renormalization-group and model calcu- 
lations In the field-theoretic calculations ||ll|,|l|, 
however, the method of dimensional regularization was 
employed and no proof was given for the unimportance 
of the cutoff dependence of fs and Cs ■ 

In this Letter we reexamine the validity of the scal- 
ing prediction Eq.(j3|) on the basis of exact results for 
the 0{n) symmetric 93* field theory at finite cutoff in 
the large-n limit. We shall show that fs depends signif- 
icantly on the cutoff procedure : a sharp cutoff A in k 
space causes non-negligible finite-size effects above and at 
Tc that violate the universal scaling form Eq.(||). Specif- 
ically, we find that Eq.(^) must be complemented as 

/.(t,L,A) = A'^-2 <i>(r^A-i) + L-'^HL/O (4) 

where the function $ has a finite geometry-dependent 
critical value <I'(0) > 0. For d > 2, the non-scaling 
term exhibits a dominant size dependence compared 
to the L^'^ scaling term. By contrast, for the (p^ lat- 
tice model with short-range interaction, we find that no 
L^^ term exists for fs{t,L,a) and that Eq.(^) is indeed 
valid except that for L :§> ^ the exponential scaling argu- 
ment of ^ must be formulated in terms of the lattice- 
dependent "exponential" correlation length [p"5|jl^ ]. 

Non-negligible cutoff and lattice effects were alre ady 
found previously for the finite-size susceptibility |]l7| , |l8| . 
Similar non-universal finite-size effects were found in sys- 
tems with subleading long-range interactions |]l9| . These 
effects, however, were restricted to the regime L » ^ 
close to the bulk limit above Tc- The new non-scaling 
finite-size effect exhibited in Eq.(^) is significantly more 
general in that it is pertinent to the entire — 
plane including the central finite-size regime ^ ^ L. In 
particular, this effect exists at Tc where 



/,(0, L, A) = K""-' $(0) + L-'' J^(0) 



(5) 
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has a non-universal leading amplitude A'^"^ *J'(0) for 
d > 2. As an important consequence, this im- 
plies the non-universality of the critical Casimir force 
^-2^rf-2 <j)^^;j^(o) in film geometry. Furthermore this 
indicates that material-dependent properties such as sub- 
leading long-range interactions (which exist in real 
systems and which do not change the universal hulk prop- 
erties) may affect the leading L dependence of fs- 

We start from the standard ip"^ continuum Hamiltonian 



V 



(6) 



where stands for {2tt) J df^k with \kj\ < A. The reg- 
ular part reads /g*^' = (ci + C2 ro)A.'^ where ci and £2 
are d dependent constants. 

In 2 < d < 4 dimensions, the finite sharp cutoff A 
causes only negligible corrections to the (leading) singu- 
lar part fl^f^ — ^"'^ of the bulk free energy den- 
sity where y^*^) is a universal amplitude . For the finite 
system, however, the A dependence turns out to be non- 
negligible. For large ^A and LA at fixed A we find that 
the finite-size free energy /i*^-* = f'^'^^ — flf'^ attains the 
form of Eq. (H) with the non-scaling contribution 



with To = roc + o-ot for the rt-component field </j(x) in the 
cubic volume V = L''' with periodic boundary conditions. 
This model requires a specification of the x dependence 
of <(5(x) at short distances. As usual, we assume that 
the Fourier amplitudes (^k = Jy d^^x {p{:x.)e'~^^^ are re- 
stricted to wave vectors k with components kj in the 
range —A < kj < A where kj — 2Tnnj/L, j = 1,2, ...,c?, 
with rrij = 0, ±1,±2,... . The question can be raised 
whether there exists a non-negligible cutoff dependence 
of the finite-size free energy density per component (di- 
vided by fc^T), 



fcube{t, Li, A) 



-n-'^L-'^ \nZ, 



cube 



(7) 



$c«6e(e"'A-l) 



6(27r; 



d-2 



dy 



dq e 



y 



L-i 

X exp [-(l + r'A-2)y] 
and with the universal scaling function 



dy 



d-l 



(12) 



(13) 



where 



Zcube{t, X, A) — 



\{d-2)/2 



exp(-ff) 



(8) 



is the dimensionless partition function. For comparison 
we shall also consider the free energy density f{t, L, d) of 
the ip"^ lattice model 



2a' 



(9) 



with a nearest-neighbor coupling J on a simple-cubic lat- 
tice with a lattice spacing d. The factor (ksT)'-^ is ab- 
sorbed in H and H. 

Essential features of the cutoff effects on the finite-size 
free energy can already be demonstrated within the sim- 
ple Gaussian model Eq.(P) with uq = which implies 



tqc — and ^ 
to 



-1/2 



Straightforward integration leads 



ro + fc')] 



A'' luTT 



L-''ln2 



(10) 



Wd{y) 
y(G) 



d/2 



E 



\rii— — OQ 



A, 



d(4 - d) 



A, 



r(3 - d/2) 
2d-2 7rd/2(d-2) 



(14) 
(15) 



We note that ^cuLi-^/O diverges in the critical limit 
L/^ ^ which is an artifact of the Gaussian model for 
cubic geometry (arising from the isolated k = term in 

AG) . 

J cube } ' 

Here we shall show that the non-scaling structure of 
Eqs.(^ and (12) is valid more generally, beyond the 
Gaussian model, for the field theory. In the limit 
n — > 00 at fixed u^n, the free energy density is [ po| 

z 2 Vouqu 



+ iL-'^5:in[A-2(x-^+k^)] (16) 



where % -'^ is determined implicitly by 



-1 



X - =ro + 4uoni-'^^(x-i+k2) ' . (17) 

k 



In the bulk limit this yields 



f(G) 



1 



A" luTT 



The bulk free energy /(, and bulk susceptibility Xfc are ob- 
tained by the replacement L~'^ J^, and the critical 
/ In [A^^ (ro + /c^)] (11) point is determined by ro = tqc = -Auqu /^.k"^ . The 
I. regular part of /{, reads /o — ci A''- — rl/{16uon). 
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For the singular part /s = / — /o of the finite-size free 
energy above and at Tc we again find the form of Eq. (H) . 
The leading non-scaling part $cM6e(C^^A^^) turns out 
to have the same form as that of the Gaussian model, 
Eq. (p^ , except that now, in the large-n limit, the bulk 
correlation length above Tc is 



(18) 



with Co = {4uonAd/[ao(4 - d)]Y'^^-^\ 

For the subleading universal scaling part we find for 

2 < rf < 4 



a 



1 fdy 



2J y 





(19) 



where P{L/^) is determined implicitly by 

oo 

4-d 



P' 



.d-2 



dy Wd{y) e-^ . (20) 



In the bulk Hmit Eqs.(|T|) and @) yield fbs = Y^''' with 

Y =id-2)Ad/[2d{A~d)]. (21) 

These results can be extended to partially confined 
L'* X oo'^^''- geometries with periodic boundary condi- 
tions in d' < d dimensions. Eq.(P remains valid where 
now the non-scaling part reads for the Gaussian model 

^±d'ir'^-') = § ^cubeir'^-') (22) 

with <i>c«fce(C^^A^^) given by Eq.(jl^). For the finite-size 
scaling function of the Gaussian model we find 



:F!,':J,iL/0^{L/0'Y(^n 



1 f dy f TT 



2 J y \y 





{d-d')/2 



X Wd'{y)e~^''/^'>"y/^^\ (23) 

For the cp'^ theory in the large-n limit the corresponding 
result is in 2 < d < 4 dimensions 



^d,d'{L/0 



A, 



2(4 - d) 



d 



where now P{L/£^) is determined implicitly by 



P 



d-2 



d-2 



47rMrf 



dy 



(d-d')/2 



X Wd'{y)e 



(25) 



We find that the non-scaling part ^d,d'{£, ^) in the 
large-n limit has the same form as the Gaussian result, 
Eq.(|2|), with i given by Eq.(|l|). 

These results have a significant consequence for the 
critical Casimir effect. In film geometry {d' = 1) the 
Casimir force is defined as [pUlOl 



Fcasrrmr = -dffamit, L, A)/dL (26) 

where the excess free energy per unit area is given by 

/|S„(i, L, A) = Lffurnit, L, A) - LMt) . (27) 

Near bulk criticality our rcsuhs, Eqs.(|2|), (|4|) and (|2£ 
yield for d' = 1 



F< 



Casimir 



(L,C,A) = i-2A''-2<I>,^i(r^A-i) 



Casimir 



(28) 



where 



Xcasrrmrix) = (d - l)Td,l{x) - X T'^^^{x) + Y x'^ (29) 

with J-'{x) — dJ-{x)/dx. Thus the critical Casimir force 
has a leading non- universal term ~ L^^, in addition to 
the subleading universal terms ^ L^'^ of previous theories 

We have also calculated the cutoff dependence of the 
finite-size specific heat C{t,L,A) above Tc. For example 
for cubic geometry and in the large-n limit we find, at 
fixed < i <C 1, the large-L-behavior 

C{t,L,k) - Cb{t)^L-^K'^-H{ik) + 0{e-^/^) (30) 

where 

vi'(eA) = So (eA)2('^-3)$,„fc,(r'A-i) (31) 

with Bo = (4 - d) (CoA)^(2-<i)/(^ _ 2f . A similar result 
holds for L'* x oo'^^'^ geometries. 

Eqs.(|), (HI) and (^) - (|l]) are the main results of 
the present paper. We have confirmed the structure of 
these results also for the ip'^ theory with finite n within 
a one-loop renormalization-group calculation at finite A 
which yields the same form of the function <I>d_(j/(C~^ A~^) 
as in Eq.(0). Thus we arrive at the general conclu- 
sion that there exists no universal finite-size scaling form 
for the leading size dependence of the free energy den- 
sity fg near criticality of confined systems with periodic 
boundary conditions, contrary to phenomenological ex- 
pectations . Clearly this raises the necessity of reex- 
amining also the universality predictions for the leading 
finite-size effects in systems with non-periodic boundary 
conditions l|-|Jll|,|2l[. 
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The sensitivity of fs{t, L,A) with respect to the cutoff 
procedure can be related to a corresponding sensitivity 
of the bulk correlation function G(x) = < ip{x)ip{0) > 
in the range |x| ^ For example, for the ip'^ con- 

tinuum Hamiltonian Eq.(^) with an isotropic sharp cut- 
off |k| < A wc find, in the large-n limit, the oscillatory 
power-law decay above Tc 



-(d+i)/2 sin [Ax - 7r(d - l)/4] 

l-F^-2A-2 



bution 



-L/C 



subleading L dependence of finite-size effects in real sys- 
tems with non-periodic boundary conditions. 

Support by DLR and by NASA under contract num- 
bers 50WM9911 and 1226553 is acknowledged. 



G'(x) = 2A'^''^{2ttxA) 

+ O [e'""^^) (32) 

for large a; = |x| ^ ^ corresponding to the existence 
of long-range spatial correlations. (A different power- 
law decay is obtained for an anisotropic sharp cutoff 
\kj\ < A.) By contrast, G(x) has an exponential decay 
for the lattice model Eq.(|^) with purely short-range inter- 
action |l5||. For the latter case we find that the scaling 
form Eq.(H) is valid provided that the second-moment 
bulk correlation length ^ in the argument of J- is re- 
placed by the lattice-dependent exponential correlation 
length |l5||l^ ]. Specifically we find for the Gaussian lat- 
tice Hamiltonian Eq.(^, at fixed t > 0, the exponential 
large-L behavior fs{t,L,a) = L~'^ ^cube{L/£,i) with 

^c«te(i/6) =rf(V2<i)''"'^/' exp(-L/Ci) (33) [7] 

where = (a/2) [arcsinh(a/2^)]~^ is the exponential 
correlation length in the direction of one of the cubic axes. 
A similar result holds in the large-n limit. Note that the 
non-universal dependence of on d is non-negligible in 
the exponent of Eq.(||) 

As pointed out by Dantchev and Rudnick the 
presence of subleading long-range interactions causes 
leading non-scaling finite-size effects on the susceptibil- 
ity similar to those for the case of a sharp cutoff in the 
regime i 3> C Clearly it would be desirable to 

extend our present analysis to the case of such interac- 
tions. On the basis of our present results we expect the 
existence of leading non-scaling finite-size effects on fs, 
for both continuum and lattice systems with subleading 
long-range interactions, not only for L 3> ^ but also in 
the central finite-size regime ^ S> i including T — Tc. 
In particular we expect a leading non-universal size de- 
pendence of the critical Casimir effect in such systems, 
different from the L^"^ scaling prediction. 

Furthermore, according to Eq.(pO[), there exists a non- 
negligible non-universal (power-law) finite-size contribu- 
tion ^ (AL)"2 also for the singular part of the finite-size 
specific heat in the region L ^ ^ where the scaling contri- 
is exponentially small for the case of pe- 



riodic boundary conditions. Corresponding non-scaling 
terms due to subleading long-range interactions are po- 
tentially important for the analysis of the leading and 
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